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1. Introduction 



It is well known that a long flexible linear polymer in a good solvent with an attractive 
short-range force between the polymer and the container wall undergoes an adsorption 
transition [jl], |2|, [|. A good model for this phenomenon is a self-avoiding random walk 
(SAW) on some lattice with an adsorbing surface (boundary). In this model each 
monomer (step) in the bulk has a Boltzmann weight x, whereas the interaction with the 
adsorbing wall is taken into account by assigning an energy e w < to each monomer 
that is found at the surface. For temperatures T higher than the critical temperature 
T a of the adsorption, the polymer is in a desorbed phase, and for T < T a it is in an 
adsorbed phase. At the adsorption transition the number M of adsorbed monomers 
scales with the total number N of monomers as M ~ N^, where (f) is the crossover 
exponent. 

The above model has been widely studied on various lattices and via a number of 
techniques (|| and references therein). For two-dimensional Euclidean lattices exact 
results were found through conformal field theory [||, [5], |] and much numerical work 
was completed using Monte Carlo |7|, §], exact enumeration || [TO], [jl], [12], |13|| , transfer 
matrix |14|, [15|], series expansion |16[ and renormalization group [17 , |18|1 techniques. 
Bouchaud and Vannimenus developed a real space renormalization group (RSRG) 
approach to study the critical adsorption of SAW on finitely ramified fractal lattices 
(which may serve as crude models for real amorphous materials). In addition to some 
other results, these authors found the exact values of the crossover exponent <fi for the 
case of two- and three-dimensional Sierpinski gaskets (SG). 

The RSRG method was applied on other fractals, in particular on other members 
of the two-dimensional SG fractal family [20], in order to find out how the crossover 
exponent changes when properties of the fractal are systematically changed. Members 
of this family can be labeled by an integer b (2 < b < oo), and for large values of b 
the underlying lattice structure becomes more and more similar to a triangular wedge, 
while its fractal dimension tends to its Euclidean value 2. Due to the extraordinary 
complexity of the underlying combinatorial problem, exact crossover exponents were 
found |2(J only for 2 < b < 9. For larger values of b the Monte Carlo renormalization 
group (MCRG) method |2]J was developed, and the limit b —>■ oo was analyzed within 
the finite-size scaling hypothesis f2~2f , with the conclusion that <fi does not tend to its 
Euclidean value <p E = 1/2 0. 

In this paper we apply the RSRG method to calculate the crossover exponent (f) for 
the previously introduced |23], nontrivial class of SAW, so-called piecewise directed 
walk (PDW) on the two-dimensional SG family of fractals. For one type (PDW1) of 
this model we were able to find exactly for the entire SG family, whereas in the other 
case (PDW2) we found exact values of <fi for 2 < b < 1000. Comparing our results with 
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the results obtained for the common SAW model on the same fractal family |2(], |21| for 
2 < b < 100 we find that for each b the values of (ft for both models are very close. This 
is not surprising since attractive wall at the fractal boundary brings anisotropy in SAW 
paths which is similar to the effects of the directedness imposed on random walks in 
PDW models. In other words, it seems plausible to accept that disallowed directions does 
not significantly contribute to the values of (ft for the SAW model (especially for large 
values of b), which means that crossover exponent (ft calculated for the PDW model(s) 
can be considered as a good approximation for the (ft values for SAW on SG fractals. 
Then, asymptotic analysis of the exact PDW results can also provide good insight into 
the crossover region b — > oo for the SAW model, which might help in resolving the 
subtle problem of the limiting behaviour of the critical properties of statistical systems 
on fractals when underlying fractal parameters approach the corresponding Euclidean 



values §§g|, U, H, |£ 

The paper is organized as follows. In section 2 we recall the definition of the PDW 
model and present the framework of the RSRG calculation of (ft for the case under 
study. In sections 3 and 4 we present our results for (ft, together with the analysis of the 
asymptotic behaviour in the limit b — > oo, in the case of PDW1 and PDW2, respectively. 
Finally, in section 5 we discuss all our findings and related results obtained by other 
authors. 



2. The models and the crossover exponent (ft 

In this section we are going to apply the RSRG method to the PDW adsorption problem 
on the SG family of fractals. We start with recalling the basic properties of well- 
known SG fractals introduced by Given and Mandelbrot |27]]. Each member of the SG 
fractal family (labeled by b) can be constructed recursively, starting with an equilateral 
triangle (generator Gi(b)) that contains b 2 smaller equilateral triangles. The subsequent 
fractal stages are constructed self-similarly, by replacing each of the b(b+ l)/2 upward- 
oriented small triangles of the initial generator by a new generator. To obtain the 
rth-stage fractal lattice G r (b), which we shall call the generator of order r, this process 
of construction has to be repeated r — 1 times, and the complete fractal is obtained in 
the limit r — > oo. The fractal dimension of the SG fractal with spatial scaling factor b 
is given by df = ln[6(6+ l)/2]/lnfe. In the case under study, it is assumed that one side 
of the fractal is an impenetrable attractive wall. 

The first type - PDW1 - of the PDW model is defined [[SJ by requiring that within 
each generator of any order of a given SG fractal the self-avoiding walker makes random 
steps in only three directions: the main one (randomly chosen), and the two neighboring 
(obtained by rotating the main direction by ±7r/3). In the case of a generator with r = 1, 
one step is related to passing through one physically present unit triangle, whereas in 
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the case of a generator with r > 1 a step should be related to the coarse-grained walk 



through the underlying fractal structure. The PDW2 type [2J] is defined in a similar 
way: within each generator of any order four directions, determined by the point of 
entering a generator, are allowed (see figure 1). 

To calculate the critical exponent by means of the RSRG approach, it is necessary 
to consider three restricted partition functions: B^ r \ B[ r \ and \ defined as weighted 
sums over all corresponding never-starting and never-ending PDWs for the given stage 
r of the iterative construction of the lattice (see figure 2). The recursive nature of the 
fractal under consideration implies the following general form of the recursion relations: 

B' = Y,B(k)B k , 



B[ = ]T B 1 (k,l,m)B k B[B™ 



k.l.m 



B' 2 = B 2 {k 1 l 1 m)B k B[B™ 



(2.1) 



k.l.m 



where we have used the prime symbol as a superscript for the (r + l)th-order functions 
and no indices for the rth-order functions. The coefficients B(k), Bi(k,l,m), and 
B 2 {k, I, m) are the numbers of ways in which the corresponding part of the PDW, within 
an (r + l)th-stage fractal structure, can be comprised of the PDWs within the fractal 
structures of the next lower order. Then, following the RSRG procedure |19|], one should 
find the symmetric fixed point B — By — B 2 — B* of the RSRG transformation (2.1), 
which corresponds to the adsorption-desorption transition at temperature T = T a , 
linearize the RSRG transformation at that point, and find the relevant eigenvalues. 
Finally, the crossover exponent is determined through the formula 

' lnA ^ (2.2) 



where A,, is 



In A, 



OB' 



OB 



B=B* 



and A^ is the larger eigenvalue of the matrix 

/ dB[ 8B[ 
A ( dBi dB 2 



8B' 2 
dB 2 



;2.3) 



(2.4) 



B=B 1 =B 2 =B* 



3. The crossover exponent </> in the case of PDW1 



It was found earlier [^3| that the first RSRG equation in (2.1) in the case of PDW1 for 
any SG fractal has the following explicit form 



b-l 



b+k 



k=0 




(3.1) 



so that \ u (|2.3|) is equal to 



where BI is the unique positive fixed point of the transformation ( |3.1|) . The equation 
( [3.1D was derived using more general partition functions X b> i(B): 



X b>i (B) = ^R bji (l)B 



where Rb,i(l) denotes number of all possible PDW1 paths with / steps in the bulk, 
starting at the origin (0, 0) and ending at the point (b, i) in figure 3, assuming that the 
weight of a step in the bulk is B. Similarly, the explicit form of the matrix A ( |2.4|) can 
be established by introducing general partition functions Y byi and Z b i x . 



Y bii (B,B u B 2 ) = ]T Q b>i (l,m,n)B l B™B 



2 ■ 



l,m,n 



Z b ^B,B x ,B 2 ) = S btl (l,™,n)B l B™B%, 

l,m,n 

where Q bt i(l,m,n), and S b> i(l,m,n) denote numbers of all possible PDW1 paths with I 
steps in the bulk, m steps at the wall, and n steps in the layer adjacent to the wall, again 
starting at the origin (0, 0) and ending at the point (b, i) in figure 3, with the additional 
assumptions that the weight of a step at the wall is B x , and in the layer adjacent to the 
wall - B 2 . One can observe that for b = 1 and b = 2 generating functions Y bt i(B, Bi, B 2 ) 
are: 

Y lfi (B, B u B 2 ) = B 2 , Y 2fi (B, B u B 2 ) = B 2 B, 
Y\,i(B, £?!, B 2 ) = B x , Y 2}l (B, B u B 2 ) = B 2 (B + B x ), 

Y 2 , 2 (B,B 1 ,B 2 ) = B 2 1 +B 2 2 B, (3.2) 

whereas for b > 2 (see figure 3(a)) the following relations are valid: 

Y bj o{B, B x , B 2 ) — B 2 B b ~ 1 , 

Y bA (B, B h B 2 ) = B 2 B h -\B x + (b - 1)B), 

i-1 

Y b+1>i {B, B u B 2 ) = Yl B'-% k (B, B u B 2 ) + BY b>i (B, B x , B 2 ), for 2 < i < b - 1, 

k=0 

6-1 

Y b+1>b {B, B u B 2 ) = B b - k Y b 4B, B x , B 2 ) + B 2 Y b , b (B, B x , B 2 ), 

k=0 

6-2 

Y b+1>b+1 {B, B x , B 2 ) = B b - k B 2 Y bik (B, B x , B 2 ) + B 1 Y b>b {B, B Xl B 2 ) . (3.3) 

k=0 
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Starting with the above equations, one can show (see [Appendix A| ) that elements An 
and A\2 of the matrix A (p.4|) are equal to: 



A n (b) 



dB[ 



dB 1 



dY h 



6,6 



dB 1 



6-2 



2E 

k=0 



1 



'6-2 N 



k+ 1 \k. 



*\2b-k-3 



(3.4) 



A 12 (b) 



8B 2 



dY, 



6.6 



dBo 



'6-2 



i 



\k=0 



k + 2\k + l \k + l 



6-3 



2k 



b-2 



(3.5) 



where * denotes that partial derivatives are calculated at the symmetric fixed point 
B = Bi = B 2 = Bl. As one can see on figure 3(b), the restricted partition function B' 2 
for the SG fractal labeled by b is 



and 



B 2 = Z b:b (B, Bi, B 2 ) 



(3.6) 



J 6+1,6+1 



(B, Bi, B 2 ) = ^ B 2 B\Xb,b-i(B) , 



i=0 



where 



X h)0 {B) =B b 



6-1 



Xb,b~i(B) = 



and 
+ 1 (b s 



k=i 



k + 1 \k I \ k — i 



\B 



b-i+k 



as was shown in ||23[| . Then, 

dZ btb 



A 21 (b) 



dBi 



6-1 

i=i 



6-1,6-l-i 



which can be put (see Appendix B ) in the following form 



6-1 



'&-r 

k 



,k 



{BD 



-l+k 



For the last element A 22 of the matrix A (|2.4|) we have 



A 22 {b) 



dZ, 



6.6 



dBo 



-S*^6-l,6-l-i(-B) 

i=0 



for i < 6, 



(3.7) 



From (3.3), with B = B X = B 2 and 1^(5, 5, 5) = X M (S), we get 



6-1 



Xb,b(B) = B 1+1 X b -i,b-i-i(B) 



i=Q 



(3.J 



(3.9) 



(3.10) 



(3.11) 
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A 22 (b) = ^X bib (B) 



so that A 2 2{b) is equal to: 

(3.12) 

Since Bl is the unique positive fixed point of the RSRG equation B' = B, where 
B' = X bib (B), from ( ^l2|) directly follows that 

A 22 (b) = 1. (3.13) 

The explicit values of fixed points Bl and eigenvalues X u were found in |23|] for 
2 < b < 1000 and for several larger values of b. Knowing Bl, it is a matter of 
straightforward calculation to find elements A n , A u , and A 21 (equations (|3.4|),(3.5), 
and ( |3.10| )), then the larger eigenvalue of the corresponding matrix A, and finally 
the crossover exponent ( |2.2| ) . In figure 4 we present the complete sequence of 0- values 
up to b = 1000, and also for b = 5000 and b = 10000. These values are compared with 
the theoretical values obtained via the asymptotic analysis for large b, which we are 
going to present now. 

For large values of b, introducing new variable k = bx the summation over k in 
equation ( |3.4|) can be replaced by integration over x. Then, by the use of the Stirling 
formula for large n, that is 

In n\ = ( n H — ^ In n — n H — In 2tt , 
V 2/ 2 

in binomial coefficients, element An (b) can be written in the form 

A ^-Tk^S^ 2b " { ' )Ax • (3 ' 14) 

where h(x) is 

h(x) — x In x + (1 — x) ln(l — x) — - (2 — x) In B% . 

Since b is considered to be very large the integral in equation ( |3.14p is treated here by 
the use of the steepest-descent method expanding the integrand around the minimum 
of h(x), which occurs at xq = 1/(1 + \ Bl), so that An(b) can be expressed as 



By similar arguments, in ref. |23|] was shown that the first equation in RSRG 
transformation (2.1) can be put in the form 



B' - — Lp- 2 ^!) / n rsifil 



for large b, where 



g(x) = xlnx + (1 — x) ln(l — x) — -(1 + x) In I? 



has its minimum at X\ = yB/(l + yB). It is trivial to show that h(xo) = g(x\) and 
h"(xo) = g"(xi), for B = B£. Then, by direct comparison of equations ( |3.15| ) and ( |3.16| ), 
one can see that 



A u {b) 



/B* b (l + JB* b ) 



3 , rMb „ ( 1 
2--(3 + V5)— + - 



(3.17) 



since B' = Bt if B = B* b , and 



R * 3-V5 35 -V51nfe 

Bb = ^r~ + 4~ — T + ° 



(3.18) 



(which was also shown in ref. ||23|| ) . 

In a quite similar way from equation (3.5) follows 



A 12 (b) = 2 l-(B* b ) b - 1 



1 



1 1 - x 



nbBT Jo x 2 



= 2^-^(3^-5)^ + 
10 b 

whereas from equation (|3.10| ) one gets 
1 r l dx 



-2b 9 ( x ) dx 



1 + JBi 



A 21 (b) 



nbBi Jo x(l - x) 



2bg(x) = 2JBt = VE-l 



V5\nb 

loir 



+o ; 



Then, the largest eigenvalue of the matrix A (|2.4j) is approximately equal to 



(3.19) 



(3.20) 



(3.21) 



for b ^> 1, where 



d = lim X^b) = - 

b— >oo Z 



3 + (41 - 8V5) 1/2 ) = 3.9037..., 



C 2 = — ((41 - 8 v / 5) 1/2 (5 -VE) + 15(3 + Vbj) = 2.2957... . 



Using asymptotic behaviour of X u , which has also been treated in |p3|| , i.e. 
A,(6) = ^^6 + ^(3v / 5-5)ln6 + --- , 



(3.22) 



together with ( p.21|) and ( |2.2|) , we finally get the following expression for the crossover 
exponent 

ln^ 



hid 
In b 



In b 



+ 



(3.23) 
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The exponent (ft given by ( |3.23| ) is represented as the dashed curve in figure 4, and one 
can see that agreement with the exact results is quite satisfactory. In the same figure, 
we have presented the lower (pi and the upper (fi u bounds 

fa = l-{d f -d a )v , <Pu = ^- (3.24) 

for the crossover exponent 0, established in a heuristic way by Bouchaud and 



Vannimenus in ||19|| . In the above bounds d s is the fractal dimension of the adsorbing 



surface, i.e. d s — 1 for our model, whereas v, the critical exponent of the end-to-end 
distance of walk in the bulk, is equal to In hj In A„. One can observe that our results 
satisfy the condition <pi < <fi < <ft u for any b, which is not the case for the SAW model 
with no constraints (see figure 3 in 



4. The crossover exponent <\> in the case of PDW2 

For the PDW2 model in the bulk it is enough to say that within each generator of 
any order of a given fractal the walker makes random steps represented only by the 
set of four vectors which corresponds to the place of his entering the generator (see 
figure 1(b)), and it is not relevant which end of the walk is starting point. This is not 
the case in the vicinity of the adsorbing wall. See, for instance, two configurations in 
figure 5, which are symmetric under change of orientation of the walk. If the wall is 
not present both of them represent PDW2 walks, and contribute the RG equation for 
the _E?( r+1 ) parameter with the same term (B^) 6 . When the wall is present these two 
configurations contribute the B 2 parameter, first with the term (B^) 3 (B2~ ) 3 , and 
second with (I?^) 3 ^^) 2 !?^, so it is obvious that the PDW2 definition in the case of 
adsorption should be completed to avoid such ambiguities. In that sense, two types of 
PDW2 walks are possible. We shall denote them by abbreviations PDW2a and PDW2b, 
and describe in the following way: 

• PDW2a model - parts of the PDW2a walks that traverse generators lying on the 
adsorbing wall are oriented in such a way that their entrance point is always on the 
wall (first example in figure 5); 

• PDW2b model - entrance point of the PDW2b walks that traverse generators lying 
on the adsorbing wall is always in the bulk (second example in figure 5). 

As in the case of the PDW1 model, it is useful to introduce analogous partition 
functions X bii (B), Y bii (B, B u B 2 ), and Z bti (B, B u B 2 ) for both PDW2 models. Functions 
Xb t i(B) were already introduced in reference P4]| , but no closed form expression for them 
were found. Nevertheless, these functions, as well as their derivatives, obey certain 
recursive relations, by the use of which one can numerically find fixed points B^ of 
the RSRG transformation B' = B (2.1) and the eigenvalues \ u {b) Q2.3|), for very large 
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values of the SG fractal parameter b [p4\ p5| . Here we shall describe the algorithm for 
numerical determining the elements of the matrix A ( |2.4|) . 

For both versions of the PDW2 model the RSRG transformation of the function B[, 
for any b, can be expressed as B[ = Y b>b (B, B\, B 2 ), so that elements An (b) and Ay^lb) 
of the matrix A ( |2.4j ) are equal to 

A n (b) = Y b ] b (B* b ), A 12 (b) = Y b 2 b (B* b ), 



where we have used short labeling 

dY h . 



1 I td\ def Ul b>i 



Y bi( B ) 



dB x 



y b 2 m 



dY h . 



dcf UI b i 



B=B 1= B 2 



dBo 



(4.1) 



(4.2) 



B=B 1 =B 2 



It is not difficult to check that for b = 1 and 6 = 2 generating functions Y by i(B, B\, B2) 
are: 

Y lfi (B, B u B 2 ) = B 2 , Y 2fi (B, B x , B 2 ) = B 2 B(1 + B), 
Y ltl (B, B u B 2 ) = B l7 Y 2>l (B, B 1: B 2 ) = B 2 (B + B 1 ), 

Y 2 , 2 (B,B 1 ,B 2 ) = B 2 1 +B 2 2 B, 

whereas for b > 2 the following relations are valid: 

6-1 



Yb+i,o(B, Bi, B 2 ) 
Yb+1,1 (B, Bi, B 2 ) 
Yb+i,i(B, Bi, B 2 ) 
Yb+i,b(B, Bi, B 2 ) 



Y,B k+1 Y btk + B b B 2 Y btb , 

k=0 

b-l 

BY bfi + ]T B k Y btk + B b - l B 2 Y b<b , 

k=l 

Y^B l -% k + Y^B k - l+1 Y Wl + B b - l B 2 Y btb , for 2 < i < b - 1, 

fc=0 k=i 

Y,B b -% k + B 2 Y b , b , 

(4.3) 



k=0 

b-l 



Yb+i,b+i(B, Bi, B 2 ) — B B 2 Y b<k + B{Y b ^ b , 

fc=0 

as one can see on figure 6(a). From these equations one obtains recursive relations for 
the derivatives Y bi {B) and Y bi (B): 



Y b \i fi {B) = B k+1 Y bk (B) , 

k=0 

Y b \i,i(B) = BY b (B) + £ B% k (B) , 

k=l 

= ]C Bi ~ k Y b ] h (B) +J2B k ~ i+1 Y b 1 i (B) , for 2 < i < b - 1, 



Yb+iAB 



k=0 
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6-1 



Yb+i,b(B) = £ B b ~ k Y b ) k {B) + BY h \{B) , 

k=0 

= E 5 6 - A+l n: fc (5) + x 6)6 (s) , (4.4) 

fc=0 

and 

Y b 2 +lfi (B) =j:B k +% 2 k (B) + B b X b , b (B), 

k=0 

y b 2 + i,i(B) = BY b0 (B) + £ 5 fc n 2 ,(5) + B b - l X h>b {B) , 

k=l 

Y b \i,{B) =^B*-%l(B)+j2B k -*+%%B) + B b -*X b , b (B), for 2 < i < b - 1, 

k=0 k=i 

Y b 2 +1>b (B) = £ B b ~ k Y b 2 k (B) + BY b b (B) + X b>b {B) , 

(4.5) 



fc=0 

i, 



6-1 



nW( 5 ) = E ^ +1 n 2 fe (£) + E ^^ fc ^(^) • 

k=0 k=0 

Then, starting from the initial values 

X lt0 (B) = Y lfi (B, B,B) = B, X 1A (B) = Y 1A (B, B,B) = B, 
Y* (B)=0, Y^(B) = 1, Y 1 2 (B) = 1, lft(fl) = 0, 

with B = B b , for any b one can iterate relations (4.3), (4.4) and (4.5), and obtain values 
of the elements A u (b) and A 12 (b) dO|) . 

One can observe that for the PDW2a model 

B' 2 = Y b>0 {B,B l ,B 2 ), (4.6) 

so that elements A 2 i(b) and A 22 {b) are equal to 

A 21 {b) = Y b ] (B* b ) , A 22 (b) = Y b 2 Q (B* b ) . (4.7) 

which can also be obtained by the use of recursion relations (4.3), (4.4) and (4.5). 
In the case of the PDW2b model 



B 2 = Z b fi(B, B 2 ) = Z bjb (B, B 1: B 2 



and 



A 2 i (6) 



dZ, 



6.0 



dB 1 



A 22 {b) 



dZ, 



6,0 



dBo 



Function Z b $(B, B\, B 2 ) satisfies the relation 



6-1 



Zb,o{B, B 1 , B 2 ) = _B 2 E B\X b ^i(B) 



(4.8) 
(4.9) 

(4.10) 



i=0 
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as can bee seen from figure 6(b). Then, from ( |4.9| ) directly follows 

A 21 (b)=J2i(B* b ) i X^ 1 , i (B* b ), (4.11) 
1=1 

together with 

A 22 (b)= b "£B*X b ^(B* b ). (4.12) 

The form of the A 22 (b) can be further simplified, since the sum on the right side of the 
last equation is equal to X b0 (Bl)/ B\ = X bb (Bl) / B b , which is direct consequence of 
relations (4.3) with B = B\ = B 2 . Finally, we obtain 

A 22 (b) = l, (4.13) 

since X b)b (B* b ) = B* b . 

In figure 7(a) we present larger eigenvalue of the matrix A (4) for SG fractals with 
b < 1000 for both PDW2 models (elements of the matrix A were calculated by the use 



of the previously found [24, |25| fixed points Bl). It seems that A^ tends to some finite 
constant when b — > oo, for both cases, which was also supported by numerical analysis 
of the obtained data. Furthermore, previous numerical analysis of PDW2 X u eigenvalues 



24] , pBJ revealed that A„ ~ b for large b, which means that crossover exponent 4> 
tends to as const/ In 6. That conclusion is in accord with figure 7(b), where calculated 
values of <fi are depicted as functions of 1/ In b. In the same figure one can also see that, 



as in the PDW1 case, Bouchaud-Vannimenus boundaries ( |3.24| ) are obeyed. 
5. Discussion 

For the sake of a better review of global behaviour of <fi as function of the fractal 
parameter b, when number of allowed directions grows from three (PDW1 model) to six 
(SAW model), in figure 8 we present our exact results for all PDW models, together 
with the exact (2 < b < 9, |2(J) and available Monte Carlo (9 < b < 100, [H) SAW 



results. For all considered b, crossover exponents are very close for all these models, 
being monotonically decreasing functions of b. More precisely, for 12 < b < 15 values of 
(f) for PDW1, PDW2b and SAW are almost the same (see the inset in figure 8), whereas 
PDW2a values slightly departure from them. For smaller values of b (b < 10) inequality 
0PDW2a < 0pdwi < 0saw < 0PDW2b holds, whereas for larger b {b > 15), one can observe 
0PDW2a < 0PDW2b < 0pdwi- As 6 becomes larger <fi values for all three PDW models 
become closer, indicating that main contribution to the value of <fi comes from walks with 
steps in three directions allowed in the PDW1 model. The fact that PDW2a values of </> 
are always smaller than PDW2b values means that the fourth direction in the PDW2a 
model (which is forbidden to the PDW1 walker) moves the walker more efficiently from 
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the adsorbing wall, than the corresponding direction in the PDW2b model (see figure 
5). Good illustration for the interplay between the effects of increasing b and changing 
the number of step directions is the relationship between PDW1 and PDW2b data: for 
smaller b values, i.e. for smaller fractal generators, there is not enough space for the 
fourth direction to move the walker away from the adsorbing wall (so 0pDW2b > 0pdwi), 
which is not the case in larger generators (0pDW2b < 0pdwi)- 

Exponent for SAW model decreases with b more rapidly than for any other model 
considered here. For smaller b exponent 0saw is close to the corresponding PDW1 and 
PDW2b values, then, as b grows, 0saw approaches PDW2a values, and starting from 
b = 80 drops even bellow them. One could expect that this trend will continue for 
b > 100, since increasing of b makes relatively more steps available in the bulk for the 
walk to spread, than on the surface. However, according to finite-size scaling prediction 
| [22|| , should tend to 1/4 when b — > oo, which means that 0saw should have minimum 
(at least one) for some b > 100, and the same conclusion holds if one accepts the 
hypothesis that 0saw tends to its Euclidean value 0£ = 1/2. There is still no firm 
evidence that either of these assumptions is correct. On the other hand, in this paper 
we have analytically proven that 0pdwi ~~ > 0, and numerically shown that 0pdw2 — * 0, 
whereas the value of for directed walk model on Euclidean lattices [28] is equal to 1/2 
(see Appendix C), as for SAW model. Since PDW walks become directed in the limit 
b — > oo, these results definitely mean that crossover exponent on fractals does not 
necessarily approach its Euclidean value and, also, give support to our expectation that 
0saw — *■ 0. Whether this expectation is correct or not, we believe that PDW models, 
being tractable for very large values of b, may serve as excellent toy models for checking 
the validity of the finite-size scaling hypothesis and the RGMC results for random walk 
models on fractals. 

At the end one can pose the question about the dependence of the crossover exponent 
on the fractal properties. There is no simple function of a limited set of the fractal 
characteristics with which can be described, however, on the basis of equation ( 3.23| ) 
and numerical analysis of the PDW2 data, we can say that behaviour of at the fractal 
to Euclidean lattice crossover is governed by the factor l/ln6 ~ 2 — df, i.e. with the 
fractal dimension. Similar behaviour was found for the PDW critical exponent v | 23|j , 
which approaches its Euclidean value v — 1 with the correction term proportional to 
l/ln&. It is interesting to note here that direct consequence of these two findings - 
relation ~ 1 — v, was also predicted by the finite-size scaling method for the SAW 



model \22\. At the same time, the lower Bouchaud-Vannimenus boundary 0/ ( |3.24j ) 
behaves as 1 — v for large b, and for all PDW models considered here exponent 
approaches 0/ when b — > oo. The same holds for SAW model on n-simplex fractals ||29|| , 
i.e. ~ (pi when n — > oo, and in the case of simple random walks on ci-dimensional 
Sierpinski gaskets is exactly equal to the lower boundary, for any d [B0]]. In contrast to 
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these cases, Bouchaud-Vannimenus boundaries are violated for the SAW on Sierpinski 
fractals [21], as well as in some other models |T0, 31, 32, 33]. Further comparative 
analysis of all these cases should be done in order to make the only known limits for the 
crossover exponent <fi more precise, which will be a matter for future study. 
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Figure captions 



Figure 1. The fractal structure G2(3) of the b = 3 SG fractal at the second stage 
of construction, with examples of (a) PDW1 and (b) PDW2 paths. Arrows represent 
the allowed step directions of the corresponding coarse-grained PDW paths. The gray 
area at the basis of G2(3) represents the adsorption wall. 



Figure 2. Schematic representation of the three restricted partition functions, for an 
rth stage fractal structure G r (b) used in the calculation of the PDW crossover exponent 
tf>. For example, Es£' represents the PDW path that starts at the G r (b) left vertex 
that lies on the adsorption wall, and exits the G r (b) at the upper vertex that lies in 
the bulk. The interior details of the G r (b) are not shown (they are manifested by the 
wiggles of the PDW paths). 



Figure 3. Illustration of the PDW1 recursion relations (3.3) and (3.7) for the 6+1 = 6 
SG fractal. The solid straight lines represent the contribution to Yq^ (a) and Zq^q (b) 
from the PDW1 paths via the legitimate paths represented by the solid curves. Arrows 
represent the possible directions of steps. 



Figure 4. Crossover exponent cf> for SG fractals with 2 < b < 1000, b — 5000, and 
b = 10000 in the case of PDW1 model (A). The dashed line is the theoretical estimate 
(3.23) of <j) obtained by the asymptotic analysis. The two solid curves represent the 
bounds (3.24) for <p. 



Appendix A. 




where 




(A.2) 



dB 



1 B=B 1= B 2 
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Figure 5. Examples for PDW2a and PDW2b configurations on the generator of the 
b=3 SG fractal. In the PDW2a case starting point A of the walk is on the adsorbing 
wall, whereas the PDW2b starting point B is in the bulk. 



Figure 6. Illustration of the PDW2 recursion relations (4.3) and (4.10). The solid 
straight lines represent the contribution to 5*6,3 (a) and ^6,o (b) from the PDW2 paths 
via the legitimate paths represented by the solid curves. Arrows represent the possible 
directions of steps. 



Figure 7. (a) The larger eigenvalue of the matrix A (2.4) for PDW2a (O) and 
PDW2b (□) model for SG fractals with 5 < b < 1000. (b) Crossover exponent for 
SG fractals with 2 < b < 1000 in the case of PDW2a (O) and PDW2b (□) model. The 
two solid lines represent the bounds (3.24) for cf) for both PDW2 models. 



The proof of equations (A.l) proceeds in four steps (i)-(iv): 

(i) Taking into account equations (3.2) one can easily verify the relations for b = 1 and 
6 = 2 

Y? fl (B) = 0, 1^(3) = 1 

Yj fi (B) = 0, Yj tl (B) = B, Y 2 ] 2 (B) = 2B, (A.3) 

which are in agreement with (A.l). 

(ii) One can also verify that for general b: 

= 0, Y b \(B) = B b -\ (A.4) 

which follows from the first two equations in (3.3) and is again in agreement with (A.l). 

(iii) Suppose that (A.l) is valid for b < m and every possible i, as well as for b = m + 1 
and % < j < m. Then, we will prove that (A.l) is valid for b = m + 1 and % — j + 1. To 
this end we use the transformation rules that follows directly from (3.3): 

C-i,o(£) = BYl (B), 

Y^ +hl (B) =BY^ (B) + BY^(B), 

y^ + i,i( B ) = E Bi ~ k Kk(B) + BY^B), for 2 < i < m, 



fc=0 



Figure 8. Crossover exponent <j> as function of 1/ ln6 for the PDW1 (A), PDW2a (O), 
PDW2b (□) and SAW (o, data taken from references [20, 21]) models on SG fractals. 
The full line is the theoretical estimate (3.23) obtained by the asymptotic analysis of 
the PDW1 data in the limit b — > oo. The inset represents magnified region 9 < b < 92. 
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m 

Y^ +1 , m+1 (B) = ]T B m+1 - k Y^ k (B) + Y m>m (B, B, B), (A.5) 

k=0 

From (A.5) we obtain 

^m+ij+i = BY^ +l j + BY^j +1 + B(l — B)Y m j . (A. 6) 

Using (A.l) for Y^ +1 •, Y^ and ■ it is a matter of straightforward algebra to show 
that 

M-E=^(-r)(v)^. (a. 7) 

(iv) If (A.l) is valid for b < m with every possible i, and for 6 = m + 1 with i < m, then 
from (A.5) one has 

y^ + i, m+ x(B) = BY' +1>m (B) + 5(1 - B)Y^JB) + X m , m (B) . (A.8) 

Again, inserting (A.l) for Y^ +1 m (B) and m (B), together with X mim (B) (see equation 
(3.8)) one finds 

by which equation (A.l) and also ( |3.4j ) has been proven for general b (since An (b) = 

The derivation of equation (3.5) can be performed in a similar way, which we are 
not going to elaborate here. 

Appendix B. 



The objective of this appendix is to derive equation Q3.10I) . Putting (3.8) into equation 
( p.9|) one immediately obtains 



A 21 (b) = (b- l)(B* b fV>-V + (B^ EE (V) ( & * - i ^ • (B ' 1} 

Introducing label f(i,k) with 

the double sum in (IB. II) can be rewritten as 



EE/(i,*) = /(M) + /(U)+ - + /(l, 6-2) 

i=l fc=i 

+ /(2,2)+ ••• +/(2,6-2) 



+ /(6-2,6-2) 
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b ~ 2 i b ~ 2 1 (h - 1\ i fb - k - ?\ 



1=1 fc=l 1=1 z + \ " / fc=l 

and, since 

which can easily be proved by induction, we finally get ( |3.1(J| ). 



Appendix C. 

For the adsorption of directed walks on homogeneous square lattice |28j it was found 
that implicit equation 

w(l-x 2 )(l + x-xw) = 1 (C.l) 

defines the value ^(w) for which (N) — > oo, when w = exp(—e w /k B T) > w c , 
where (N) is the average lenght of the walk (x is statistical weight of a step, whereas 
e w < is an energy assigned to each step at the surface). Values x = x c = v2 — 1 
and w = w c = 1 + l/y/2 correspond to the adsorption transition. From ( |C.1| ) 
straightforwardly follows that 

w - w c ~ (xoo(w) - x c ) 1/2 , (C.2) 



for w — > and x — > i m (w). It was also shown in |£8| that 



(N) ~ , \ and (M) ' d%M 



Xaaiw) — X Xnoiw) — X 6.W 



where (M) is the average number of steps at the surface. Then, by the use of (|C.2|) , one 
easily obtains 

(M) ~ (iV) 1 / 2 , 

which means that crossover exponent is equal to 1/2 for directed walk adsorption on 
2d Euclidean lattices. 
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